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Abstract 

Direction relations between extended spatial objects are important 
commonsense knowledge. Recently, Goyal and Egenhofer proposed a for- 
mal model, known as Cardinal Direction Calculus (CDC), for representing 
direction relations between connected plane regions. CDC is perhaps the 
most expressive qualitative calculus for directional information, and has 
attracted increasing interest from areas such as artificial intelligence, ge- 
ographical information science, and image retrieval. Given a network of 
CDC constraints, the consistency problem is deciding if the network is re- 
alizable by connected regions in the real plane. This paper provides a cubic 
algorithm for checking consistency of basic CDC constraint networks, and 
proves that reasoning with CDC is in general an NP-Complete problem. 
For a consistent network of basic CDC constraints, our algorithm also 
returns a ' canonical' solution in cubic time. This cubic algorithm is also 
adapted to cope with cardinal directions between possibly disconnected 
regions, in which case currently the best algorithm is of time complexity 
0(n 5 ). 



1 Introduction 

Representing and reasoning with spatial information is of particular importance 
in areas such as artificial intelligence (AI), geographical information systems 

*This paper is an extended version of Xiaotong Zhang, Weiming Liu, Sanjiang Li, Ming- 
sheng Ying: Reasoning with Cardinal Directions: An Efficient Algorithm. AAAI 2008: 387- 
392. 
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(GISs), robotics, computer vision, image retrieval, natural language process- 
ing, etc. While the numerical quantitative approach prevails in robotics and 
computer vision, it is widely acknowledged in AI and GIS that the qualitative 
approach is more attractive (see e.g. [4]). 

A predominant part of spatial information is represented by relations be- 
tween spatial objects. In general, spatial relations are classified into three 
categories: topological, directional, and metric (e.g. size, distance, shape, 
etc.). The RCC8 constraint language [26] is the principal topological formal- 
ism in AI, and has been extensively investigated by many researchers (see e.g. 
^9\^7\^U^2lM\^B^2\^2\- When restricted to simple plane regions, RCC8 
is equivalent to the 9-Intersection Model (9IM) [6], which is a very influential 
relation model in GIS. 

Unlike topological relations, there are several competitive formal models for 
direction relations [7, 8 , 2 . Most of these models approximate a spatial object by 
a point (e.g. its centroid) or a box. This is too crude in real- world applications 
such as describing directional information between two countries, say, Portugal 
and Spain. Recently, Goyal and Egenhofer [TTJ [10] proposed a relation model, 
known as cardinal direction calculus (CDC), for representing direction relations 
between connected plane regions. In CDC the reference object is approximated 
by a box, while the primary object is unaltered. This means, the exact geometry 
of the primary object could be used in the representation of the direction. This 
calculus has 218 basic relations, which is quite large when compared with RCC8 
and Allen's Interval Algebra pQ. Due to its expressiveness, CDC has attracted 
increasing interest from areas such as AI [32j [33j [23] , GIS [12] , database [31] , 
and image retrieval [T4] , 

One basic criterion for evaluating a formal spatial relation model is the 
proper balance between its representation expressivity and reasoning complex- 
ity. While reasoning complexity of the point-based and the box-based model of 
direction relations have been investigated in depth (see [20] and [2]), there are 
few works discussing the complexity of reasoning with CDC. 

One central reasoning problem with CDC (and any other qualitative cal- 
culus) is the consistency (or satisfaction) problem. Other reasoning problems 
such as deriving new knowledge from the given information, updating the given 
knowledge, or finding a minimal representation can be easily transformed into 
the consistency problem [4]. In particular, given a network of CDC constraints 

J\f = {vi^ijVj^j^i (each Sij is a CDC relation) (1) 

over n spatial variables , v n , the consistency problem is deciding if the 

network M is realizable by a set of n connected regions in the real plane. The 
consistency problem over CDC is an open problem. Before this work, we do 
not know if there are efficient algorithms deciding if a set of CDC constraints 
are realizable. Even worse, we do not know if this is a decidable problem. 
Furthermore, given a satisfiable set of CDC constraints, how to construct a 
realization in the real plane? 

This paper is devoted to solving these problems. We first show by examples 
that local consistency in particular path-consistency is insufficient for deciding 
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the consistency of basic CDC constraints. Then we devise a cubic algorithm 
for checking if a network of basic CDC constraints is consistent. In case the 
network is consistent, this algorithm also generates a solution that is canonical 
in a sense. Moreover, we also show that deciding the consistency of an arbitrary 
network of CDC constraints is an NP-Complete problem. This implies that 
reasoning with CDC is decidable. 

Some restricted versions of the consistency problem have been discussed in 
the literature. Cicerone and di Felice [3] discussed the pairwise consistency 
problem, which decides when a pair of basic CDC relations (<S, 5') is consistent. 
Skiadopoulos and Koubarakis [32 investigated the weak composition problem 
[5| IT8] of CDC, which is closely related to the consistency problem of basic CDC 
networks involving only three variables. 

The CDC algebra is defined over connected regions. A variant of CDC was 
proposed in [33], where cardinal directions between possibly disconnected re- 
gions are defined in the same way. This calculus, termed CDC^ in this paper, 
contains 511 basic relations. An 0(n 5 ) algorithm was proposed in [33 for check- 
ing the consistency of basic constraints in CDC^, but the consistency problem 
over CDC is still open. Recently, Navarrete et al. [23] tried to adapt the ap- 
proach used in [33] to cope with connected regions, but their approach turns 
out to be incorrect. 

The remainder of this paper proceeds as follows. Section 2 recalls basic no- 
tations in qualitative spatial/temporal reasoning and introduces the well-known 
Interval Algebra (IA) [1 and Rectangle Algebra (RA) — the two-dimensional 
counterpart of IA. We introduce the CDC algebra in Section 3, where the connec- 
tions between CDC and RA relations are established in a natural way. Section 4 
introduces the notion of canonical solution of a consistent basic CDC network. 
Section 5 first proposes an intuitive 0(n 4 ) algorithm for consistency checking of 
basic networks and then improves it to 0(n 3 ). We apply our main algorithm to 
the pairwise consistency problem and the weak composition problem over CDC 
in Section 6. In Section 7 we adapt the main algorithm for connected regions 
to solve consistency checking in two variants of CDC. Conclusions are given in 
the last section. 

2 Qualitative Calculi: Basic Notions and Exam- 
ples 

Since Allen's Interval Algebra, the study of qualitative calculi or relation models 
has been a central topic in qualitative spatial and temporal reasoning. This 
section introduces basic notions and important examples of qualitative calculi. 

2.1 Basic Notions 

Let Dbea universe of temporal or spatial or spatial-temporal entities. We use 
small Greek symbols for representing relations on D. For a relation a on D and 
two elements x, y in D, we write (x, y) G a or xay to indicate that (x, y) is an 
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instance of a. For two relations a, /3 on D, we define the complement of a, the 
intersection, and the union of a and (3 as follows. 



We write Rel(D) for the set of binary relations on D. Clearly, the 6-tuple 
(Rel(D); — , n,U,0,D x D) is a Boolean algebra, where and D x D are, 
respectively, the empty relation and the universal relation on D. 

A finite set B of nonempty relations on D is jointly exhaustive and pairwise 
disjoint (JEPD) if any two entities in D are related by one and only one relation 
in B. We write (B) for the subalgebra of Rel(Z}) generated by B, i.e. the 
smallest subalgebra of the Boolean algebra Rel(D) which contains B. Clearly, 
relations in B are atoms in the Boolean algebra (B). We call (B) a qualitative 
calculus on D, and call relations in B basic relations of the calculus. A similar 
definition was given by Ligozat and Renz [21] , where B was required to be closed 
under converse and contain idjj — the identity relation on D. 

For two relations a,^on D, the converse of a and the composition of a and 
(3 are defined as usual. 

a ~ = {(2/5 x) e D x D : (x,y) e a} 
a o (3 = {(x, y) G D x D : (3z G D) [(x, z) G a and (z, y) G (3}}. 

The composition of two relations a,f3 in (B) is not necessarily in (B), i.e. ao (3 
cannot be represented as the union of some relations in B. We say a qualitative 
calculus (B) is closed under composition if the composition of any two relations 
in (B) is still a relation in (23). In general, for a, (3 G (B), the weak composition 
[5j [18j [28] of a and /?, written as a o w (3, is defined to be the smallest relation 
in (B) which contains a o (3. Clearly, a qualitative calculus is closed under 
composition if and only if the weak composition operation is the same as the 
composition operation. 

An important reasoning problem in a qualitative calculus (B) is the consis- 
tency (or satisfaction) problem. Let A be a subset of (B). A constraint over A 
has the form (xjy) with 7 G A. For a set of variables V = {vi}f =1 , and a set of 
constraints M involving variables in V, we say AT is a constraint network if for 
each pair there exists a unique constraint (vijVj) in M. A network M is 
said to be over A if each constraint in M is over A. In particular, we say a net- 
work is a basic network if it is over B. A constraint network AT = {^i7ij^j}^j=i 
is consistent (or satisfiable) if there is an instantiation {ai}^ =1 in D such that 
(di, dj) G jij holds for all 1 < i, j < n. In this case, we call {ai}^ =1 a solution of 
J\f . The consistency problem over A is the decision problem of the consistency 
of constraint networks over A. 



aD (3 
aU/3 



—a 



{(x,y) G D x D : (x,y) a} 

{(x, y) G D x D : (x,y) G a and (x, ?/) G /3} 

{(x, 2/) G L> x D : (£, 2/) G a or (x, y) G (3}. 
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2.2 Interval Algebra 



Interval Algebra (IA) [1 is a qualitative calculus defined on the set of (closed) 
intervals in the real line. I A is generated by a set Bi nt of 13 JEPD relations 
between intervals (see Table HI. 



Relation 


Symbol 


Converse 


Meaning 


before 


P 


Pi 


x~ 


< x + < y~ 


< y + 


meets 


m 


mi 


x~ 


< x + = y~ 


< y + 


overlaps 





oi 


x~ 


<y~ < x + 


< y + 


starts 


s 


si 


x~ 


= y~ < x + 


< y + 


during 


d 


di 


y~ 


< x~ < x + 


< y + 


finishes 


f 


fi 


y~ 


< x~ < x + 


= y + 


equals 


eq 


eq 


x~ 


= y~ < x + 


= y + 



Table 1: Basic IA relations and their converse, where x = [x , y = [y , y + ] 
are two intervals. 

IA is closed under converse and composition. This means IA is a relation 
algebra in the sense of Tarski [35] . The computational complexity of reasoning 
with IA has been extensively investigated by researchers in artificial intelligence 
(see H5] and references therein). In particular, Allen [1 introduced the 
important notion of path-consistency for networks of IA constraints and proved 
that path-consistency suffices to decide the consistency of basic IA networks. 
Nebel and Biirckert [25] found a maximal tractable subclass of IA, which is 
known as ORD-Horn and usually denoted by TL. For a network of IA constraints 
over 7Y, path-consistency also suffices to decide consistency. 

The definitions of basic IA relations as given in Table [l] concern only the 
ordering of the endpoints of intervals. This suggests that different solutions of 
the same basic IA network respect the same ordering. In particular, we could 
choose intervals that have integer endpoints. 

Definition 1 (canonical set of intervals [22]). Suppose m = { [m^ , m^]}^ =1 is 
a set of intervals. Set E(m) to be the set of endpoints of intervals in m. We say 
m is a canonical set of intervals iff E(m) = [0, M] D Z, where M is the largest 
number in E(m). A solution of a basic I A network is called a canonical interval 
solution if it is a canonical set of intervals. 

Clearly, if m = {[m~ , m^]}^ =1 is a canonical set of intervals, then each m~ 
(mf) is an integer between and 2n — 1. Moreover, set M to be the largest 
number in E(m). Then M < 2n and for any < m < M there exists i such 
that m~ = m or = m. The following theorem shows that each consistent 
basic IA network has a unique canonical solution. 

Theorem 1. Suppose M = {viXijVj}fj =1 is a basic IA network. If AT is con- 
sistent, then it has a unique canonical interval solution. 

Proof. Suppose {h}™^ is a solution of A/*, where li = [^ r ,Z^]. Write ao < 
ai < ••• < a n * for the ordering of {^~,^}£=i- Define f,g : {!,-•• ,n} — > 
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{0, 1, • • • ,n*} as f(i) = k if l~ = otk and g(i) = k if if = a^. Let hi = 
[f(i),d(i)]- Because only the ordering of endpoints of intervals matters in a 
solution, it is easy to see that () = {hi}f =1 is a solution of AT. Set E(§) to be 
the set of endpoints of all intervals hi. It is easy to see that n* = max 
Moreover, there exists i such that m = f(i) or m = g(i) for each < m < n*. 
In other words, E(§) = [0,n*] D Z. Therefore, \) is a canonical interval solution 
of M. Such a solution is clearly unique. □ 

2.3 Rectangle Algebra 

Rectangle Algebra (RA) [13j |2] is a qualitative calculus defined on the set of all 
rectangles in the plane, where we assume a predefined orthogonal basis in the 
plane and only consider rectangles or boxes two sides of which parallel to the 
axes of the orthogonal basis. 



I y (r l )o\I y (r 2 ) 



rip <g> oir2 



Figure 1: Illustration of Rectangle Relation 

For a rectangle r, write I x (r) and I y (r) as, resp., the x- and ^/-projection of 
r. The basic rectangle relation between two rectangles ri,r 2 is defined by the 
basic I A relation between I x (r\) and I x (r 2 ) and that between I y (r{) and I y {r2) 
(see Fig. [I]). More precisely, if (I x (ri), I x (^2)) G ® and (I y (ri), Iyfo)) G /?, then 
we write a ® (3 for the basic rectangle relation between r± and r 2 . In other 
words, for any basic IA relations a, (3, 

(ri,r 2 ) Gft®^ (4(ri),4(r 2 )) G a and (/ y (ri), 7 y (r 2 )) G /?. (2) 

Write B rec for the set of these rectangle relations, i.e. 

B rec = {a®/3:a,l3eB in t} (3) 

RA is defined as the qualitative calculus generated by B rec on the set of rect- 
angles. 

Proposition 1. A basic RA network AT = {vi^ij <S> f3ijVj}fj =1 (otij,(3ij G Bi nt ) 
is consistent iff its component IA networks M x = {^a^i^-}™ - =1 and M y = 
{vi(3ijVj}™j =1 are consistent. 



6 



Proof. This naturally follows from the definition of basic RA relations. □ 

The definition of canonical solution can be easily extended to RA. 

Definition 2 (canonical rectangle solution). For a consistent basic RA network 
J\f = {vidij <g> f3ijVj}™ j=1 (aij,Pij G B int ), a set of rectangles {m i }^ =1 is a 
canonical rectangle solution of M iff its x- and ^/-projections, {I x (^i)}i=i an d 
{I y (mi)}i=ii are canonical interval solutions of M x = {viaijVj}fj =1 and M y = 
{vi0ijVj}? j=1 , respectively. 

3 Cardinal Direction Calculus 




— • ^ l*i=^ H 2 =b 8 ; b^ = b 9 



Figure 2: A bounded connected region b and its 9-tiles 

In this section we first introduce Cardinal Direction Calculus (CDC) of Goyal 
and Egenhofer [11] and then establish its connection with Rectangle Algebra. 
In particular, we will associate a basic RA network Af r with each basic CDC 
network A/*, such that M is consistent only if M r is consistent. More profound 
connection will be established in the following section. 

3.1 Direction Relation Matrix 

CDC is a qualitative calculus defined for extended two-dimensional objects in 
the real plane. 

Definition 3 (plane region and connected region). A subset a of the real plane 
is called a region if a is a nonempty regular closed subset, i.e. if a = a°, where x° 
and x are the (topological) interior, and respectively, the (topological) closure 
of a subset x of the real plane. A region a is said to be connected if it has a 
connected interior a° . 

For a bounded set b in the real plane, set 

x~(b) = inf{x : (x,y) G 6}, x + (b) = sup{x : (x,y) G 6}, (4) 
y~(b) = w£{y : (x,y) G 6}, y~(b) = sup{y : (x,y) G b}. (5) 
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We write 



4(6) = [x-(b),x+(b)}, I y (b) = [y-(b),y + (b)}. 



(6) 



Set 



M(b) = I x (b) x I y (b). 



(7) 



We call M(b) the minimum bounding rectangle (mbr) of 6, and call I x (b) and 
I y (b), respectively, the x- and ^/-projection of b. Clearly, Ai(b) is the smallest 
rectangle containing b. 

Remark 1. The RA relations can be extended from rectangles to arbitrary 
bounded regions via their minimum bounding rectangles. The extended RA 
relation between two bounded regions a and b is defined to be the RA relation 
between M(a) and M(b). 

By extending the four edges of M(b), we partition the plane into nine tiles, 
denoted as NW b ,N h ,NE b , W b ,O b ,E b , SW b , S b ,SE b (see Fig. |2] (left)). For ease 
of representation, we also write in sequence &n, &12, • " ' ^33 or b 1 ,b 2 , • • • , b 9 for 
these tiles (see Fig. [2] (right)). Note that each tile is a region, and the intersection 
of two tiles is either empty or of dimension lower than two. 

Since the partition only concerns the mbr of 6, the following lemma is clear. 

Lemma 1. For two bounded regions b,c, if M(b) = M(c), then b % = c % for 
i = 1,2, ••• ,9. In particular, if we set m = Ai{b), then b l = m l for each 
z = l,2,-.- ,9. 

The following notion of direction relation matrix was first proposed by Goyal 
and Egenhofer [TT] for representing the cardinal direction between connected 
plane regions. 

Definition 4 (direction relation matrix). Suppose a, b are two bounded con- 
nected plane regions. Take b as the reference object, and a as the primary object. 
The direction of a to b is encoded in a 3 x 3 Boolean matrix dir(a, b) = [dij]ij =1 , 
where = 1 if and only if a°Dbij 7^ 0, where a° is the interior of a (see Fig.[3|. 
In this case, we call dir(a,6) a direction relation matrix, or a valid matrix. 

By Lemma [l] it is clear that the direction of a to b is the same as that of a 
to M(b). 

Lemma 2. For connected regions a^b^c, if Ad(b) = A4(c), then dir(a,b) = 
dir(a,c). In particular, dir(a,b) = dir(a,M(b)). 

This lemma shows that the shape of the reference object is irrelevant: what 
matters is its mbr. 

In what follows we make no distinction between a valid 3x3 Boolean matrix 
and the direction relation it represents. It is easy to see that not all matrices 
are valid. That a matrix is valid is closely related to the following notion of 
4-connectedness. 

Definition 5 (4-connected Boolean matrix). Amxn Boolean matrix [<iij]i<i< m ,i<j< 
is said to be 4-connected if for any two nonzero cells pq and st, there are a series 
of k + 1 nonzero matrix places PoQo^PiQir ' ' ^PkQk such that 
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• pq = PoQo and st = p k q k \ and 

• Pi+iQi+i is a 4-neighbor of p^, i.e. \p i+ i - pi\ + - = 1, for each 
z = (),••• 

In other words, two cells in a Boolean matrix are 4- neighbors if they are 
horizontally or vertically adjacent, and a Boolean matrix is 4-connected if any 
two nonzero cells are connected by a series of 4-neighbors. 

Proposition 2. fWjj A 3 x 3 Boolean matrix is a direction relation matrix if 
and only if it is nonzero and 4-connected. 

Goyal and Egenhofer identified all together 218 valid matrices. For a pair of 
bounded connected region a, 6, it is clear that (a, b) determines a unique direc- 
tion relation matrix, viz. dir(a, b). Write Bdir f° r the set of cardinal directions 
represented by these matrices. This means, Bdir is a JEPD set of relations. 
Cardinal Direction Calculus (CDC) is defined to be the qualitative calculus 
generated by Bdir over the set of bounded connected regions. 




/ \ / 1 1 \ / 1 1 \ 

dir(6,o) = dir(6,a / ) = 1 1 , dir(o,6)=[ 1 1 ,dir(a / ,6)= 1 

\ 1 1 / \ / \ / 

Figure 3: Illustrations of basic CDC relations 

In order to describe the relative position of two connected regions a, 6, know- 
ing the direction of b to a is not enough. Fig. [3] shows an example, where the 
direction of b to a is the same as that of b to a' but the direction of a to b is dif- 
ferent from that of a' to 6, i.e. dir(6, a) = dir(6, a') but dir(a, b) ^ dir(a', b). This 
is drastically different from IA and many other well-known qualitative calculi. 
For example, the basic I A relation of a to b is uniquely determined by that of b 
to a. Section 6.1 will investigate this pairwise consistency problem in detail. 

In the remainder of this section, we establish the connection between CDC 
and RA relations. 

3.2 Direction Relation Vector 

We begin with the one-dimensional counterpart of CDC relations. 

Definition 6 (direction relation vector). Suppose / = [x~ , x + ) and J = [y~ , y + ) 
are two intervals. Interval J partitions the real line into three parts L\ = 
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(— oo, L2 = (?/ _ ,2/ + ), and L3 = [?/ + ,+oo). The direction of J to J is 
encoded in a Boolean vector dir(7, J) = (di, <i 2 , ^3), where ^ = if and only if 
n Li = 0. In this case, we call (^1,^2,^3) a direction relation vector. 

Clearly, a Boolean vector t = (£1,^2, £3) is a direction relation vector if and 
only if there exist two intervals /, J such that t = dir(7, J). The following lemma 
gives a characterization of direction relation vectors. 

Lemma 3. A Boolean vector t = (^1,^2,^3) is a direction relation vector if and 
only ift^ (0,0,0) and t ^ (1,0,1). 

Interestingly, each direction relation vector actually represents an IA rela- 
tion. 

Lemma 4. For two intervals /, J, suppose t = (£1, £2^3) is the direction relation 
vector of I to J. Then we have 

• t= (1,0,0) iff IpJ orlmJ; 

• t = (0, 1, 0) iff IsJ or Id J or If J or IeqJ; 

• t = (0, 0, 1) iff I pi J or ImiJ; 

• t = (1,1,0) iff IoJ orlfiJ; 

• t = (0,1,1) iff Jo/J or I Si J; 

• t = (1,1,1) iff IdiJ. 

Proof Take £ = (1,0,0) as example. For two intervals / = and J — 

[y~, dir(J, J) = (1,0,0) if and only if x~ < y~ , (ar,z+) n (^,?/ + ) = 0, 
and x + < y + hold. This is equivalent to saying that x + < y~ , which is possible 
if and only if IpJ or ImJ. □ 

In what follows, we call an I A relation a vector I A relation if it is the I A 
relation represented by a direction relation vector. We make no difference be- 
tween a direction relation vector and the vector IA relation it represents. By 
the above lemma, we know there are six vector I A relations, viz. 

p u m, s u d u f u eq, pi u mi, u fi, oi u si, di 

Note that a vector IA relation is in general non-basic, but a pair of vector IA 
relations are more precise. For example, from dir(7, J) = (0,1,0), we are not 
sure whether Is J, or Id J, or if J, or /eq J hold. Assuming dir( J, /) is also given, 
then it is easy to see that the IA relation between /, J is definite, i.e. a basic 
IA relation. 

The following lemma summaries the correspondence between pairs of direc- 
tion relation vectors and I A relations. 

Lemma 5. For a pair of direction relation vectors (s,t) and two intervals /, J, 
we have s = c//f(/, J) and t = dir( J, I) if and only if (/, J) is an instance of one 
basic IA relation in the cell specified by (s,t) in Table^ 

Lemma [5] shows that all basic IA relations except 'meets' and 'before' (and 
their converses) can be represented as pairs of direction relation vectors. 
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s \ t 


(1,0,0) 


(0,1,0) 


(0,0,1) 


(1,1,0) 


(0,1,1) 


(1,1,1) 


(1,0,0) 








p,m 











(0,1,0) 





eq 





f 


s 


d 


(0,0,1) 


pi, mi 

















(1,1,0) 





fi 














(0,1,1) 





si 





oi 








(1,1,1) 





di 















Table 2: Pairs of vector I A relations. 



3.3 Projective IA Networks 

The connection between CDC and RA relations is established via the notion of 
projective interval relations. 

Definition 7 (projective interval relation). For a basic CDC relation <$, the 
x-projective interval relation of S is defined as 

i*(5) = {(I x (a),I x (b)):(a,b)€6}, (8) 

where I x (a) and I x (b) are the x-projective intervals of M(a) and A4(6), respec- 
tively. 

The following lemma proves that each t x (S) is indeed a vector IA relation. 

Lemma 6. Suppose S = [dij]^j =1 is a basic CDC relation. Then the x-projective 
interval relation l x (5) is the IA relation associated to the vector (^1,^2,^3), i-e. 
(/, J) G t x (S) iff dir(I,J) = (d 1 ,d 2 ,d 3 ), where dj is if T% =1 dij = and 1 
otherwise. 

Proof. See Appendix A. □ 

Therefore, each x-projective relation t x (S) is an IA relation. Immediately, 
we have 

Lemma 7. For a pair of basic CDC relation (5,7), if (a, b) is an instance of 
{v\ 5v2, ^27^1}, then (I x (a), I x (b)) is an instance of l x (S) and the converse of 
**( 7 ), i.e. (4(a), I x (b)) e t x (S) n **( 7 )~ 

Proof. From (a, b) E 5 and (6, a) £ 7 we have (I x (a), I x (b)) G ^(5) and (I x (b), I x (a)) G 
^(7). Therefore, (I x (a),I x (b)) £ l x (5) C\ l x . □ 

As a consequence, we have 

Lemma 8. ^4 kmc CDC constraint network M = {^Aj^j}™j=i ^ consis- 
tent only if the IA constraint network {vitfjVj}" j =1 is consistent, where t x j = 

Proof. Suppose {a i }^ 1 is a solution to M. Then {I x (a>i)}i=i 18 a solution to 

{n'fjrAfj^. '"" ' □ 
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Note that ifj is empty or a basic IA relation or the non-basic IA relation 
p U m or its converse. Set 

B* nt = {o, s, d, f, eq, fi, di, si, oi} u {p u m, pi u mi} (9) 

For each nonempty IA relation i in B* nt , define 

{p, if^ = pum; 
pi, if l = pi U mi; (10) 
otherwise. 

We call Tthe meet-free refinement of i. Clearly, T= t\(mU mi) for each i in 

u mt' 

Lemma 9. An I A network M = {vi^ijVj}i'j=i over B* nt is consistent if and 
only if M = {viiijVj}™j = i is satisfiable. 

Proof See Appendix B. □ 

All the above notions and results also apply to the ^-direction. 

Definition 8 (projective IA networks). For a basic CDC network M = {viSijVj} 7 j L j z 
recall is an IA relation defined in Lemma [8| We write M x and M y , resp., for 
the basic IA networks {vipfjVj}™j = i and {vipfjVj}™j =1 , where 

p% = ^.\(mumi) (11) 
Pi = ^\(mumi). (12) 

We call M x and M y , respectively, the x- and ^/-projective I A networks of M. 
We write M r for the basic RA network {viPijVj}™j=n where pij = pfj p\-. 

As a corollary of Lemma [8] and Lemma [9j we know 

Theorem 2. A basic CDC network of constraints M = {viSijVj}™j =1 is consis- 
tent only if the projective IA networks M x = {^iPfj v j}^j=i andM y = {viPijVj}™j =1 
are consistent. 

The above theorem shows that if a basic CDC network is consistent, then 
its projected IA networks are also consistent. By Prop. [I] this implies that the 
associated basic RA network M r is also consistent. 

Example 1. Fig. [^] specifies a basic CDC network M — {vi$ijVj}ij =1 , and its 
associated RA network M r = {^iPij^j}ij=i, where pij — pfj <S> p\y For each 
pair of i ^ j, a solution of {viS i} and a solution of {vipijVj} are also 

illustrated. 

In the next section, we prove that a consistent CDC network M has a solution 
{ a i}r=i sucn that {M(ai)}^ =1 is a solution to M r , the associated basic RA 
network of M. 
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Figure 4: A basic CDC network and its associated basic RA network 



3.4 Local Consistency Is Insufficient 

As we have mentioned in §2.2, path-consistency is sufficient to decide the consis- 
tency of a basic IA network. This property is shared by several other qualitative 
calculi, including RCC8 [24, 29 and the point-based cardinal direction calculus 
of Ligozat [20]. This property, however, does not hold for CDC. 

We first note that a basic network is path-consistent if and only if every 
subnetwork involving three variables is consistent. The following example shows 
that there is a basic CDC network that is path-consistent but still inconsistent. 
Similar conclusion has been obtained in [33] for the case of possibly disconnected 
regions. 

Example 2. Let a\, a 2 , as, be the squares illustrated in Fig. [5| Consider 
the basic CDC network AT = {^ij^j}ij=i, where 

• Sij = dir(ai, dj) for 1 < ij < 4, 





• So* 



1 




for 1 < i < 4, and 





1 








1 




i 





; 




























1 1 

• S 61 = I 1 ) ,5 52 = S 53 = 


Fig. ^and Fig.^show that every subnetwork of N which involves four variables 
is satisfiable. This implies in particular that J\f is path- consistent. 

M is, however, inconsistent. We prove this by contradiction. Suppose 
{bi}i =1 is a solution to N ' . Then M(bi)(l < i < 4) will be related in a configu- 
ration similar to that of ai in Fig. [5[ By the assumption of 5^ for 1 < i < A, 
we know Ai(bi) is contained in Ad(b^) for each i = 1,2,3,4. By the definition 
of S51, we further conclude that M.(b\) contains the lower left corner of M.(b$). 
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Similarly, Ad(b^) contains the lower right corner of M.(b$). Meanwhile, accord- 
ing to 5$i we have b§ D M,(bi) = for 1 < i < 4. By the configuration of 
&i , &2 , &3 , ^4 (c./- [5p, £Ws means no point 0/65 can 6e feeZow £fte line y = yo, 
where yo = sup{n : (x, n) G M.{bi)}. This is impossible since M.(b\) C M.(b§). 
Therefore, JV is inconsistent. 

In fact, for any positive integer /c, the consistency of all subnetworks of AT 
with k variables does not guarantee the consistency of JV. A counter example 
can be constructed analogously (see Fig. [7| where k squares instead of four are 
used). This shows that local consistency is insufficient for consistency even for 
basic CDC networks. To determine if a CDC network is consistent, we need to 
consider the constraints in a network as a whole. 

4 Canonical Solution 

Suppose JV = {vi5ijVj}™j=i is a consistent basic CDC network, where each Sij 
is in Bdir- We show JV has a canonical solution in a sense similar to that for I A 
network. We begin with an arbitrary solution a = {ai}™ =1 of A/", then transform 
regions in a step by step into regions in the digital plane 1? without changing 
the cardinal direction relations between any two of these regions. 

4.1 Regular Solution 

Suppose a = {a^}-^ is a set of n connected regions. This subsection shows how 
to regularize a. Illustrations are given in Fig. [8] 

Write mi = [x^,xf] x [y~[ ,yf] for the mbr of a^. Set 

e~ = min{x~ : 1 < i < n}, (13) 

e + = max{x+ : 1 < i < n}, (14) 

/~ = min{n~ : 1 < i < n}, (15) 

/ + = maxj^ : 1 < i < n}. (16) 

Let 

5(a) = [c-c+]x [/"/+]. (17) 

















a 2 


a, 




a 4 



Figure 5: A solution to the subnetwork over {i>i, v 2 , ^3, ^4} of JV in Example [2] 
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Figure 6: Solutions to other subnetworks of J\f in Example exam:path- 
consistencyNOTcons which involve four variables 









b 3 






b 2 











b k 



Figure 7: A configuration of k squares 




Figure 8: Illustration of regularization 
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Extending edges of each rectangle ra$ until meeting the boundary of 5(a), we 
partition 5(a) into small cells. Suppose 

a < a\ < • • • < a Ux (18) 

is the ordering of real numbers in {x~ ,xf : 1 < i < n}, and 

00 < Pi < ■ ■ ■ < Pn y (19) 

is the ordering of real numbers in {y~ , yf : 1 < i < n}. Denote 

Cij = [otiiOti+i] x (0 < i < n x ,0 < j < rty) (20) 

and write 

C(a) = {cij : < i < n x , < j < n y } (21) 
for the set of these cells. For each z, let 

a\ = \J{c e C(a) :cno^ 0}. (22) 

Definition 9 (frame, cell set, regularization). For a set of connected regions 
a = {ai}f =1: we call 5(a) the frame of a, and call C(a) the cell set of a, and 
call a r = {a[}^ =1 the regularization of a, where 5(a), (7(a), and a\ are defined 
in Equations |17|21|22l respectively. 

The regularization a r has the same frame and the same cell sets with a. 

Lemma 10. For a set of connected regions a = {di}2 = i, we have ai C a\ 
for each i, and 5(a) = 5(a r ) and (7(a) = C(a r ), where a r = {a[}^ =1 is the 
regularization of a. 

Proof Straightforward. □ 
Example [l] (continued) 

For the constraint network specified in Fig. HJ Fig. [8] illustrates how to trans- 
form a solution {ai, <22, as} (leftmost of Fig. to a regular solution {a^, a^, ag} 
(rightmost of Fig. |8|. The frame and cell set are also illustrated. 

If a is a solution to a basic CDC network A/", then so is its regularization. 

Lemma 11. Suppose a = {ai}™ =1 is a solution to a basic CDC network N '. 
Then a r = {a[}^ =1 is also a solution to N '. 

Proof It is clear that each a\ is connected, and M(d[) = M.(di). It is straight- 
forward to show that dir(a^,aj) = dir(a[,ap for any z, j. □ 

Definition 10 (regular solution). A solution a = {ai}^ =1 of a basic CDC net- 
work M is called regular if a is the same as its regularization a r = {a[}^ =1 . 

It is easy to see that the regularization of a solution is regular. This is 
because a r and a have the same frame and the same cell set. By Eq. 22, it is 



straightforward to see that a solution a = {a^}^ =1 of a basic CDC network J\f is 
regular if and only if each is the union of all cells in (7(a) that has nonempty 
intersection with the interior of a,-. 
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4.2 Meet-Free Solution 



We next show each consistent basic CDC network J\f has a solution that is 
meet-free in the following sense. 

Definition 11 (meet-free solution). A solution a = of M is meet- free 

if for any z, j, I x (di) does not meet I x (aj), and I y (cLi) does not meet I y (aj). 



Suppose a = is a solution of J\f and a r = {a[}^ =1 is its regularization 

(see Fig. 11). Suppose rrii = [x^,xf] x [y~r ,yf] and rrij = [xj,x~j~] x [yj ,yf\ 
meet at x direction, i.e. xf = xj . Recall cyq < ■ ■ ■ < a Ux is the ordering of 



: 1 < i < n}. We call = x~\ an x-meet point. Clearly, k > and 



x~ < CKfc-i < x^ = xj = ak < ot-k+i < We next show how to delete this 
meet point by transforming a r into another regular solution. 

Write a* = (afc + afc + i)/2. The line x = a* divides each cell (0 < Z < n y ) 
into two equal parts, written in order and c^. For each 1 < s < n and each 
< / < n yi if Cki C a£ but c^-i^i 52 a s then delete cj^ from a£. The remaining 
part of each a£, written as 6 S , is still connected, and it is straightforward to show 
that b = {b s }™ =1 is also a regular solution of J\f. Such a modification introduces 
no new meet points. Continuing this process for at most n times, we will have 
a solution that has no x-meet points. The same modification can be applied to 
7/-meet points. In this way we obtain a meet-free solution. An example is shown 
in Fig. [9] It is easy to see that the meet-free solution has the same frame as a. 
The two solutions, however, have different cell sets. 











a{ 






— a{ 


V 

























l 2 a" 



f 2 a * 



Figure 9: Illustration of meet-freeing 

In conclusion, we have the following lemma: 

Lemma 12. Each consistent basic CDC network has a regular solution that is 
meet-free. 

So we can safely assume that a is a meet-free solution that is regular, i.e. 
a r = a. 
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4.3 Digital Solution and Canonical Solution 

Suppose a is a meet-free solution of M that is regular. We next transform a into 
a solution in the digital plane Z 2 . 

Definition 12 (pixel, digital region, digital solution). A pixel is a rectangle 
Pij = + 1] x + 1], where z, j are integers. A region a is digital if a is 
composed of pixels, i.e. D a° =^ iff p^- C a. A solution a = {ai}™ =1 of a 
basic CDC network is digital if each is a digital region. 



Pl2 



■ i ■ i 

i i < i 



■a 




■ i i ■ 



1 £ 3 4 5 6 



Figure 10: A pixel and a digital region a 



As in the paragraph immediately above Lemma 11 , we write mi = [x. 



[yiiVt] for tne mbr of a i- Recall 5(a) and C(a) = {cij}o<i<n x ,o<j<n y are, 
respectively, the frame and the cell set of a. Since a is regular, each is 
composed of a subset of cells in C(a). 

We next show how to transform the solution a = {a s }™ =1 into a digital 
solution b = {b s }™ =1 . For 1 < s < n, define a subset b s of S(a) = [0, n x ] x [0, n y ) 
as follows: 

A pixel p^ is contained in b s iff the cell Cij in C(a) is contained in a s . 
That is, 

b s = [J{Pij : (kj Q Q>s, < i < n>x, < j < %} (23) 



Clearly, 6 S is a connected digital region (See Fig. 11 for illustration). By def- 
inition we have dir(a^aj) = dir(^,6j) for any Therefore, the assignment 
b = {6 s }g =1 is also a solution of J\f. We observe that b has several good prop- 
erties. 

Lemma 13. Suppose a = {a s }™ =1 is a regular meet-free solution of a basic 
CDC network AT. The digital solution b = {b s }™ =1 constructed as above is also 
regular and meet-free. Moreover, {I x (b s )}™ =1 ctnd {I y (b s )}™ = i are the canonical 
solutions of the projective I A networks J\f x and J\f y , respectively. 

Proof. That b is regular and meet-free follows directly from the properties of 
a. We show {I x (b s )}™ =1 and {I y (b s )}™ =1 are canonical sets of intervals. Take 
the x-direction as example. Because b — {bi}f =1 is a solution to A/*, we know 
(I x (bi), I x {bj)) is an instance of t x (Sij) D L x (5ji)~. Since b is meet-free, I x {bi) 
cannot meet I x (bj) for any two z, j. This implies that {I x (bi)}f =1 is a solution 
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Figure 11: Transform a regular solution into a digital one 



to the basic IA network Af x = {^iPf 7 ^j}™ 7= i, where is defined by Eq. fTTl By 

icalsel 



the choice of (Eq. 18), we can easily show {/ x (^)}iLi is a canonicalset of 



intervals. □ 

We call the solution b in the above lemma a canonical solution. More pre- 
cisely, we have 

Definition 13 (canonical solution). A solution b = of AT is said to be 

canonical if it is digital, regular, and meet-free, and {/ :E (6i)}^ =1 and {I y (bi)}f =1 
are canonical sets of intervals. 

As a corollary of Lemma [l3j we have 

Theorem 3. Each consistent basic CDC network has a canonical solution. 

Canonical solutions are not necessarily unique. For a consistent basic CDC 
network JV = {vi5ijVj}™j =1 , recall we write M x = {viPfjVj}™-! and M y = 
{viPijVj}fj =1 for the two projective basic IA networks (see Dfn.|8|. By Thm.[2j 
we know both J\f x and J\f y are consistent. Suppose X = {Ii}2 =1 and J = { Ji}2=i 
are their canonical interval solutions as given by Thm. u\ where Ii = [x^,xf], 
Ji = [y^,y^~]- Write n x = max{x^}- l =1 < 2n, n y = max{^ + }- l =1 < 2n. Set 
T = [0,nj x [0,n y ] and mi = Ii x J^. Clearly, C T. By Thm. [3| M has a 
canonical solution because of its consistency. 

Lemma 14. If a — {a^}^ =1 is a canonical solution of JV, then each ai is a 
connected digital region whose mbr is rrii = Ii x Ji, where {Ii} f =1: { Ji}f =1 are, 
respectively, the canonical interval solutions of M x and J\f y . 

As a corollary, we have 

Corollary 1. Suppose b = {bi}™ =1 and b' = {^}™ = i are two canonical solutions 
of AT. Then I x (b % ) = I x (b[), I y (b z ) = I y (bfi, and M(b % ) = M^). 

This shows that two canonical solutions have the same mbrs. Moreover, 
it can be proved that the union of two canonical solutions is also a canonical 
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solution. Is there a maximal canonical solution? Where a canonical solution 
& = { a i}?=i of a basic CDC network M is said to be maximal if for any other 
canonical solution a' = {a-}™ =1 of AT we have a\ C di for z = 1, • • • , n. 

In the next subsection, we give a method for constructing the maximal canon- 
ical solution. 

4.4 Maximal Canonical Solution 

For a consistent basic CDC network A/*, suppose a = {ai}™ =1 is a canonical 
solution of Af. By Lemma [14| and the definition of canonical solution, we know 
{I x (ai)}i=i and {I y (Q>i)}i=i are the canonical interval solutions of, respectively, 
Af x and J\f y . Set li = I x {a>i), J% = I y {a>i), and mi = Ii x for each z. Because 
of the constraints in A/*, not all pixels in can appear in 

Note a basic CDC relation S = [d st \l t= i can be rewritten as a 9-tuple S = 
{ ( ^ ) )%=i by setting d^ = d st , where <p = 3(s — 1) + 1. For each variable Vj and 
each 1 < (j) < 9, we write ra^ for the 0-th tile associated to the rectangle rnj 
(see Fig. [i] (right)). For each constraint 5ij = (dfj)^ =1 of Vi to Uj, by {a^aj) 
satisfies Sij and .M(aj) = m^, we know d^- = if and only if a° D = 0. 

Therefore, each pixel contained in is disallowed in ai if d^- =0. For each z, 
we write 

= {Pst ^ : there exist j 7^ z, (j) such that C 777^ and d^- = 0} (24) 

That is, Di contains all pixels in that are disallowed for violating some 
constraint in Af. Set 

bi = [j{Pst £ m» : p 3t £Di} (25) 

Lemma 15. If a = {a^}™ =1 zs a canonical solution of Af , then M.{ai) — AA{bi) 
and a° P\p s t — for any p st G i.e. a$ C bi, where bi is defined in Eq. \25^ 

We note that bi may be disconnected. As a connected region, must be 
contained in a unique connected component of bi. We observe that two con- 
nected digital regions share at least one pixel in common if they have the same 
mbr. 

Lemma 16. Let a and b be two connected digital regions. If a° D b° — 0, then 
M(a)^M(b). 

Proof. Without loss of generality, suppose A4(a) = [0, n\] x [0, 77,2], where 77i, 772 
are positive integers. Suppose a° D b° = and A4(a) — A4(b). Since A4(a) = 
[0, Tii] x [0, 77*2], a contains a pixel po s and a pixel p ni -i,t for some s, t G {0, 1, • • • , 
n2 — l}. Because a is a connected digital region, there is a path a (i.e. a sequence 
of pixels in a) that connects po s to Pm-i,t- Clearly, a separates M(a) into at 
least two disjoint components, each of which is contained in a rectangle that 
is smaller than A4(a). Since a° D b° = 0, we know 6 is contained in one such 
component. Therefore, A4(b) is smaller than A4(a). A contradiction. □ 
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By the above lemma, we know each bi has at most one connected component 
Ci such that A4(ci) = m^. Moreover, if Af is consistent, then each bi has a unique 
component q with M(ci) = M(bi) = rn h . 

Lemma 17. Suppose Af = {viSijVj}fj =1 is a consistent basic CDC network. 
Then each bi has a unique connected component Ci such that M(ci) = M(bi) = 

Vfli. 

Proof. Suppose {a^}™ =1 is a canonical solution. We have is connected and 



Ai(ai) = rrii. By Lemma 15 , we know C Let c$ be the unique component 
of bi which contains a$. Then by m % = AAiaf) C A4(ci) C A4(bi) = rn t we know 
M(ci) = rrii. □ 

Moreover, {q}™ =1 is also a solution of Af ' . 

Lemma 18. Let C{ be the unique connected component ofbi such that M(ci) = 
rrii. Then c = {ci}^ =1 is a solution of Af . Moreover, ifa= {a$}£=i ^ s a canonical 
solution of N ' , then ai C a for each i. 

Proof. That is contained in q is clear. We need only show dir(Q,c J ) = Sij 
also holds for each pair of z, j. This is equivalent to proving for all 1 < <\> < 9 
that c° fl = if and only if d^- = 0. Recall ai C q C ^ C = .M(a^). 

If d^ =0, by q C ^, for each pixel p st G we know c° fl^ = 0, hence 

c ? ° H raf = 0. On the other hand, if df, = 1, then c° H mf D a ? ° H raf ^ 0. This 
shows c is a solution of Af. □ 

It is easy to see that each region q in c is a connected digital region. More- 
over, c is a regular and meet-free solution of Af such that I x (ci) = Ix(p>i) an d 



Iy( c i) — Iy( a i) f° r each i. By Dfn. 13 , we know c is also a canonical solution of 
Af. By ai C a and the arbitrariness of a, we know c is the maximal canonical 
solution of A/". 

Theorem 4. Suppose Af = {vib'ijVj}™ - =1 zs a consistent basic CDC network. 



Let Ci be the unique connected component of bi of Eq. 25 such that A4 (q 



Then c = {ci}f =1 is the maximal canonical solution of Af. 

It is worth stressing that each m^, ^, q are independent of the choice of 
canonical solution a. They only depend on the specific constraint network Af. 

An example is given in Fig. [l2] (left) to illustrate the procedures. Note that 
each bi obtained in this example happens to be connected. 

Example [l] (continued) 

For the network specified in Fig. [4j we have 

< x^ < X2 < x% < xf < £3" (26) 



The canonical interval solutions of Af x and Af y are illustrated in Fig. 12 (left). 
Suppose {a, 6, c} is a solution of the network described in Fig. [4] Note that the 
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(2,2)-entry of 5 12 is 0, which is possible only if a° D Ai(b) — 0. This excludes 
pixel P12 from a (see Fig. 12). Fig. [l2| (right) illustrates the maximal canonical 
solution {a, 6, c} of A/". 




Figure 12: Canonical interval solutions and maximal canonical solution 



5 A Consistency Checking Algorithm 

In this section, we describe our algorithm for checking consistency of basic CDC 
networks. As in the last section, we assume J\f = {vAjVj}^j =1 is a basic CDC 
network. To examine the consistency of A/*, we first compute the two projective 
I A networks AT X and Af y . If any of these two networks is inconsistent, then 
J\f is inconsistent, either. Assume both M x and AT y are consistent. We then 
compute their canonical interval solutions, and therefore construct a frame T 
and a rectangle ra$ for each i. As in the last section, we continue to compute 
the digital region bi and try to find the connected component q such that 
A4(ci) = rrii. If such Ci does not exist for some z, then M is inconsistent. 
Otherwise, set c = {q}^ =1 . If c is a solution, then M is consistent. Otherwise, 
as implied by Thm. [4j Af must be inconsistent. 
Fig. [13] gives the flowchart of the algorithm. 

5.1 An 0(n 4 ) Consistency Checking Algorithm 

In this section, we give a detailed description of our consistency checking algo- 
rithm. 



Step 1. Projective IA Networks 



For any basic CDC constraints {xi5i2X2, ^2^21^1}, the x- and ^/-projective IA 
relations pf 2 and p\ 2 (Eq.s 11 and 12) can be computed in constant time. So the 
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1 




Figure 13: Flowchart of the main algorithm 



projective IA networks M x and M y can be constructed in 0(n 2 ) time. If M x or 
My is inconsistent, which can be checked in cubic time, then M is inconsistent. 

Step 2. Canonical Interval Solutions 

Suppose M x and M y are consistent. Their canonical solutions and 
{Ji}f =1 can be constructed in cubic time. 

Write Ii = [x^,x^~] and Ji = [y^,y^]- By definition of canonical interval 
solution, we know x^ , x\ , , yf are integers between and 2n — 1. Write 
n x = max{x^} and n y = max{^ + }. Set T = [0,nj x [0,n y ] as the frame and 
define m, = [x^,xf] x [?/~,^ + ] for each i. 

Step 3. Excluding Impossible Pixels 

As before, we write D{ for the set of pixels in ra$ that are disallowed for i.e. 

Di = {Pst ^ ™>i : there exist j ^ i, such that p st C and d^- =0} 
Equivalently, p st G ^ if and only if p s£ C mi D |J{ m j : ^ = 7^ 0- Set 

Clearly, b° D is empty for any j, with d^- = 0. Note that each bi can 
be computed in 0(n 2 ) time from Di. To compute Di, an intuitive method is 
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checking for each pixel p in ra$, and each tile rrij of any reference object rrij, 
whether p is in rrij. Note that there are at most 0(n 2 ) pixels contained in rrii 
and at most 0{n) different rrij . This requires cubic time for a fixed z, and hence 
0(n 4 ) time in total. Later, we will show this can be improved to 0(n 3 ). 



Step 4. Connected Components 

We further compute connected components of bi for each i. Note that bi has 
at most one component whose mbr is rrii (Lemma 16). If no such component 
exists for some z, then Af is inconsistent. Otherwise, set q for the connected 
component of bi such that M(ci) = rrii. 

Applying a general Breadth-First Search algorithm, we can find all connected 
components of bi and determine if their mbrs are rrii m 0{n 2 ) time. So, this 
step need only 0(n 3 ) time. 



Step 5. Checking A Possible Solution 

The last step is then checking if {q}™ =1 is a solution of Af. Note that if the 
answer is yes, then {q}™ =1 is the maximal canonical solution of Af . 

For each pair q and Cj, we should check if dir(Q,Cj) = Sij. In other words, 
we should check for each 1 < <fi < 9 whether the following equation holds 

c° H = 4» 4 = ( 28 ) 

When dfj =0, by the definition of bi, we have b° flm^ = 0. Since q is contained 
in bi, we know c° D rrij = 0. The condition is always satisfied. 

So to determine if dir(Q,Cj) = S^, we need only check whether c° D is 
nonempty for each <fi with d^. = 1. 

Recall Ci is a connected component of bi and A4(ci) = ra$. If m° D rrij = 0, 

then c° fl raf = 0. The condition is violated and hence Af is inconsistent. Note 

that each rrii has been computed and each rrij can be computed in constant 

time. Therefore, whether m° D rrij is empty can be checked in constant time for 
any z, j,^. 

Suppose rrii° D 7^ 0. Write ra^ D = x [y~ ,y~^~]. To check if 

c° D rrij is nonempty, we need show there is a pixel p which is contained in both 
rrii H and q. We need not check this for all pixels in rrii H ra^ . Instead, we 
need only check whether p/^ C q for all boundary pixels p^i of D rrij with 
(fc,Z) G Hi U i?2 7 where 

= {(fc,Z) : fc e {x~ ,x + - 1} and y~ <l < y + }, 
H 2 = {(M) : ^~ < < x+ and Z g {y~,y + - 1}}. 

We justify the above statement as follows. If m^Dra^ = rrii, then by A4 (q) = 
we know there exists a boundary pixel which is contained in q. Otherwise, 
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rrii DmJ is a rectangle strictly contained in m^, which means q contains a pixel 
p out of in rrii D raff. If q contains no boundary pixel of rrii D m^, then, because 

it is connected, q contains no pixel of rrii H at all. 

Since i^i Ui72 contains O(n) pixels and checking if a pixel is contained in q 
needs constant time, dir(Q,c J ) = 5^ can be checked in O(n) time. 

In conclusion, we can determine in cubic time whether {q}™ =1 is a solution of 
A/*. Now, since only Step 3 needs at most 0(n 4 ) time, the algorithm determines 
the consistency of a basic CDC network in 0(n 4 ) time. 

5.2 A Cubic Improvement 

In this subsection, we improve the main algorithm to cubic. This is achieved 



by an 0(n 2 ) improvement for computing each Di (see Eq. 24) in Step 3. As a 
consequence, all bi can be computed in cubic time. 

Suppose T = [0, n J x [0,n y ]. Then any digital region a C T can be repre- 
sented as an n x x n y Boolean matrix B(a) as follows: 

b <«>[m - { I: — . <»> 

Example 3. The left of Fig. illustrates a digital region a contained in 



T = [0,4] x [0,4], and the right of Fig. 14 shows the Boolean matrix B(a) 
that represents a. 
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Pu 
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Pn 








Poi 










Poo 








POA 



B(a) 



/0 1 1 1 0\ 
1 
1 

\0 0/ 



Figure 14: A digital region a and the Boolean matrix B(a) 
It is worthy noting that we adopt a different way to address the elements of 



B(a). For example, the lower left corner of B(a) in Fig. 14 is addressed as the 
(0,0)-element, i.e. 5(a) [0,0], instead of the (5,1) -element. 

For each z, recall Di is the set of pixels in rrii that are disallowed for Vi 



(cf. Eq. 24 or Step 3 of the main algorithm). For simplicity, we write Pi for the 
matrix that represents the digital region which consists of pixels in Di. Similarly, 
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we write Bj for the matrix that represents the rectangle m^, and write for 

the matrix that represents Tflra^. In case TDm^ contains no pixel, we assume 
it is represented by the zero matrix. 

Boolean operations A and V can be defined on Boolean matrices in a natural 
way. By definition of Di (see Eq. [24]) , we have 



Pi = BiA \/{Bf : 4 = and j + i}. (30) 
Adding up all Bj with df- = and j ^ z, we obtain an integer matrix Qi, 

i. e. 

Qi = Y,i B f ■ 4 = and 3 + i}- (31) 
Easy to see, the following equation holds: 

P .\ k n = { if QA k A > and Bi[k,l] = 1; 
l[ ' J \ 0, otherwise. ^ ' 

Note that T D is a (may be degenerative) rectangle. We next show that 
this property plays an important role in the improvement. To this end, we 
introduce the following operations on matrices. 

Definition 14 (cumulative matrix and difference matrix). For a matrix TV, we 
define its cumulative matrix as 

k 

acc(JV)[M] =5^JV[t,Z]. (33) 
t=o 

and its difference matrix as 

dW.q = {^;J'_^ [Jb _ 1>q> otwt. ( 34 ) 

The cumulative matrix can be computed column by column. We first add 
the first column to the second one, and then add the updated second column to 
the third, etc. In this way, cumulative matrix acc(TV) can be computed linearly 
in the number of elements of N. 

It is easy to verify that N = acc(diff(7V)). Since the difference operation is 
additive, i.e. diff(7Vi + N 2 ) = diff(TVi) + diff(7V 2 ), we have 

Qi = acc(diff(Q0) 

= acc(diff(5^{B/ : 4 = and J + *})) 
= acc(^{diff(s/) : 4 = and J + *}) 



We next show each Qi can be computed in quadratic time. To this end, we 
need the following lemma, which asserts that the number of non-zero elements 
in diff(Kf) is of order 0(n). 
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Lemma 19. The number of non-zero elements in diff(Bf) is less than 4n. 



Proof. If Bj is the zero matrix, then diff(I^) is also the zero matrix. Otherwise, 
the non-zero elements in Bj compose a rectangle, i.e., there exist < x\ < X2 < 
n x and < y\ < 2/2 < n y s.t. 



Bf[k,l] 
Easy to prove, 

diff(2#)[fc,Z] = 



r i, if 

1 0, ot 



x\ < k < X2 and y\ < I < y<i\ 
otherwise. 



1, if k = x\ and y\<l <yi\ 

-1, if k = X2 + 1 and yi < I < y2\ 
0, otherwise. 



So if ^2 < n x — 1, there are (2/2 — 2/1 + 1) 'l's and '-l's in diff(£>J); otherwise 
#2 = — 1, there are (?/2 — 2/1 + 1) 'l's and none '-l's, with other elements being 
zeros. So there are at most 2 x (2/2 — 2/1 + 1) < 2 x n y < 4n non-zero elements 
indiff(^). □ 

Since 1 < j < n and 1 < (j) < 9, there are at most 9n Bj need to consider. 
By Lemma 
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we know ^{diff(I^ ) : dfj = 0} can be computed in quadratic 
time. As it is~a matrix with n x x n y elements, computing its cumulative matrix, 
viz. Qi, also needs quadratic time. 

By Eq.[32j we know P^, the matrix representation of Z^, can be computed in 
quadratic time. As a consequence, each hi can be computed in 0(n 2 ), instead of 
0(n 3 ), time. In this way, we improve Step 3 of the main algorithm from 0(n 4 ) 
to cubic time. 

Example [l] (continued) 

In our running example, for mi, there are four different such that Bj is not 



the zero matrix and 



n 2 



0, viz. ?77| 



m2, ^25 m 2' m 3- We have 
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\ 1 1 1 ) 



diff(pt) 
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V i o o 
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-1 
-1 







Summing up all corresponding diff(Pj) and calculating its cumulative ma- 
trix, we obtain Q\. 
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By Eq. 32 


we know Pi[fc, I] 


= min{Si[fe,Z],Qi[fe,Z]}. 


Therefore 



Pi = min{ 



/ 1 1 1 \ 
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110 
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V J 



Finally, we get b\ by Eq. [25j shown in Fig. [l2j which is clearly connected. 
This shows c\ =b\. 



5.3 Beyond Basic CDC Constraints 

In the above two subsections, we have shown that the consistency of a basic 
network of CDC constraints can be determined in cubic time. Using backtrack- 
ing method, we immediately know that the consistency satisfaction problem of 
CDC is an NP problem. 

Lemma 20. The consistency satisfaction problem of CDC is an NP problem. 

Proof. Let C = {viCijVj}^^ be a set of CDC constraints. To determine if C is 
consistent, we need only branch each non-basic constraint c^-, and then call our 
cubic algorithm to solve the basic network of CDC constraints. □ 

This problem is also NP-hard. 

Lemma 21. The consistency satisfaction problem of CDC is NP-hard. 

Proof. We prove this by reducing a known NP-hard problem to the consistency 
satisfaction problem of CDC. Let A be the JEPD set of I A relations 

{p u m u pi u mi, o u s u d u f u eq u fi u si u di u oi}. 

It can be proved that reasoning with A is already NP-hard. Actually, the NP- 
hardness of A is guaranteed by the work of Krokhin, Jeavons, and Jonsson 
[15] . We need only to show A is not contained in any of the eighteen maximal 
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tractable subclasses of IA (see [I5j Table III]). The verification is straightforward 
and omitted here. 

Reasoning with IA relations in A can be easily reduced to reasoning with 
CDC. For a set of I A constraints 

AT = {vi\ijVj}™ j=1 , (X i:j e A), 

it is easy to see that J\f is satisfiable if and only if the set of CDC constraints 

A/** = {vi(pijVj}t j=1 

is consistent, where cpij is the disjunction of the basic CDC relations in 
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if = p U m U pi U mi, and is the disjunction of the basic CDC relations in 

/ \ / \ / \ / \ 
{ 1 0,1 1 0,0 1 1,1 1 1 }, 
\ / \ / \ / \ / 

otherwise. 

This shows that reasoning with CDC is at least as hard as reasoning with 
A. Therefore, reasoning with CDC is also NP-hard. □ 

As a corollary, we know 

Theorem 5. Reasoning with CDC is an NP-Complete problem. 

6 The Pairwise Consistency Problem and The 
Weak Composition Problem 

Our main algorithm can be applied to solve the two special subproblems of 
reasoning with CDC, i.e. the pairwise consistency problem and the weak com- 
position problem. These two subproblems have been considered in [3] and [32] . 
respectively. In this section, we will compare their results with ours. 

6.1 The Pairwise Consistency Problem 

Given that you know the relation of a to 6, what about that of b to a? Math- 
ematically speaking, this is the converse problem. The relation of b to a is the 
converse of that of a to 6, and vice versa. Suppose B is a set of JEPD relations 
on D. The qualitative calculus (B) is not necessarily closed under converse. 
This means, a~ could be a relation outside of (B) despite a G (B). 
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As for CDC, we have shown in Fig. [3] that a basic CDC relation may have 
more than one 'converses,' where for two basic relations a,/?, we say a is a 
converse of (3 (in CDC) if {viav 2 ^v 2 Pvi} is consistent. 

The pairwise consistency problem in CDC is the problem of deciding if 
{v\5i 2 v 2l ^2^21^1} is satisfiable for a pair of basic CDC relations Si 2 and S 2 i. 
This problem has been discussed by Cicerone and di Felice [3] , where they iden- 
tified 2004 consistent pairs of basic CDC relations. 

We next apply our main algorithm to solve the pairwise consistency problem. 
The first step computes pf 2 = l x (S 12 ) H l x (5 21 )~ and p\ 2 = L y (S 12 ) H t y (5 21 )~. If 
pf 2 or p\ 2 is empty, then the program stops and returns 'inconsistent.' Other- 
wise, we go to Step 2. We construct the canonical solutions {Ji, 1 2} and { Ji, J 2 } 
to p x 2 and p y 2: respectively. Set = Ii x Ji for i = 1,2. Write Ii = [x~ ,xf) 
and Ji = [y~,y^]- Denote n x = maxjx^x^}, n y = max^,^}. Clearly, 
1 < n X: n y < 3. Set T = [0,nj x [Q,n y ]. Then T C [0,3] x [0,3]. Continuing 
as described in the main algorithm, we will determine if {vi5i 2 v 2 ,v 2 5 2 iVi} is 
consistent, and find the maximal canonical solution in case it is consistent. 

A specialized algorithm is implemented. We obtain in total 757 consistent 
pairs of basic CDC relations. Among these consistent pairs (£, 5'), one S may 
correspond to multiple 5' ', i.e. 6 may have multiple converses. In fact, among 
the 218 basic CDC relations, 119 have unique converse, 68 have two converses, 
6 have four converses, 20 have eight converses, 4 have thirty converses, and one 
(viz. the one such that dij = 1 iff i = j = 2) has 198 converses. 

Our result is unexpectedly different from that of [3 , where Cicerone and di 
Felice obtained 2004 consistent pairs of basic CDC relations. A careful exami- 
nation, however, shows that a similar but different model was used in [3]. 

When defining the direction relation matrix dir(a,6) = [dij]ij =1 of a to b 
(see Dfn. [4|, we require = 1 if and only if a° D bij 7^ 0. While in [3], the 
(r/)-entry d^ is 1 if and only if a has nonempty intersection with 6^-, i.e. d^ = 1 
iff aDbij 7^ 0. We call these two definitions the interior-based and, respectively, 
the closure-based direction relation matrix. 

In the following, we argue that the interior-based definition is more coherent 
than the closure-based one. 

First, though the original definition of Direction Relation Matrix (DRM) [11 
does not mention that d^ = 1 iff a has a common interior point with the tile b^ , 
Goyal and Egenhofer defined in [12] the detailed direction relation matrix of a 
to b as a numerical matrix dir*(a, b) = [d*j]fj=i, where d*j is interpreted as the 
ratio of the area of aDbij and a. This is a natural extension of the interior-based 
direction relation matrix. From dir*(a,6) = [^*j]f J= i, we can obtain the coarse 
direction relation matrix dir(a, b) = [dij]\ j =1 by setting d^ — 1 iff d*j > 0. 

Second, the qualitative calculus introduced by the interior-based direction 
relation matrix is more desirable. For example, it is a natural requirement 
that the identity relation is contained in a unique basic CDC relation. For the 
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interior-based definition, we have 



/ \ 
dir(a,a) = 1 J . 

\ / 

for any connected region a. The following example, however, shows that this is 
not the case for the closure-based direction relation matrix. 

Example 4. Take a as the square [—1,1] x [— 1, 1], and take a' as the unit disk 
centered at (0, 0). Then M(a') = Ai(a) = a, but c//r(a, a) ^ dir(a' , a') if we take 
the closure-based definition. 

/ 1 1 1 \ / 1 

dir(a,a)= 111, dir(a',a') = 1 1 1 
\ 1 1 1 / \ 1 



Moreover, the interior-based definition is also consistent with the one adopted 
in [32j [33]. 

6.2 The Weak Composition Problem 

The notion of weak composition plays a very important role in qualitative spatial 
and temporal reasoning [TJ [5j [18j [28] . For two basic CDC relations a,/?, the 
weak composition a o w (3 of a and (3 is defined to be the smallest relation in the 
CDC algebra which contains the composition a o (3. Since CDC is a Boolean 
algebra, a o w f3 is the union of all basic CDC relations it contains. For a basic 
CDC relation 7, it is easy to prove that 

7C ft oJO7n(fto^)/0. (35) 

Note that 7 D (a o f3) is nonempty iff the following set of basic CDC constraints 

C = {v!av 2 , v 2 (3v 3 , V!jv 3 } (36) 

is consistent. We note that C is not a complete network. The constraint of, say, 
V2 to i>i, is not specified. According to the previous subsection, a may have 
multiple converses. To apply our main algorithm, we need to extend C to a 
complete network: 

C* = {^1^2, v 2 a'v 1: v 2 f3v 3 , v 3 (3'v 2 , ^17^3, ^37^1} (37) 

We then call our main algorithm to determine if the above completed network 
is consistent. If the answer is 'yes' for some a' ', f3' ',7' ', then 7 is contained in 
the weak composition of a and (3. Note that we need only to apply the main 
algorithm to those a 7 , /?', 7' such that (a, a 7 ), (/?, (3') : (7, 7') are consistent pairs. 

The above algorithm was implemented, and its codes are available from the 
authors. 
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The same problem has been considered in [32], where Skiadopoulos and 
Koubarakis gave an algorithm to compute the weak composition. The main 
idea is to compute the weak composition progressively. 

Definition 15 (single-tile, multi-tile, component). A basic CDC relation is 
called single-tile if its matrix has only one nonzero entry, and called multi-tile 
otherwise. We say a single-tile relation [sij]^j =1 is a component of a multi-tile 
relation [dij]f }J - = i, if Sij < dij for 1 < i,j < 3. 

The weak composition of two single-tile relations has been computed by 
Goyal [10 . Upon this, Theorem 1 of [32] establishes a rule for computing the 
weak composition of a single-tile relation and a basic relation. The correctness 
of this theorem is confirmed by our algorithm. Furthermore, Theorem 2 of [32] 
then gives a rule to compute the weak composition of two multi-tile relations. 

Definition 16 (decomposability). Let a, (3 are two basic CDC relations. Sup- 
pose ai, Q2? • * * i a k are the component single-tile relations of a. We say a basic 
relation 7 is decomposable with respect to (a, (3), if 7 is the join of k direction 
relation matrices 7 S , where each 7 S is a basic relation contained in a s o w (3. 

The following rule was used in [32 to compute the weak composition between 
a and any other basic CDC relation (3. 

\32\ Theorem 2] A basic relation 7 is contained in a o w (3 if and only if 7 is 
decomposable with respect to (a, (3). 

We next give two examples to show that this rule is incorrect in some cases. 
Example 5. For basic relations a,/?, 7 in Eq. 



(38) 



we assert that 7 is not decomposable with respect to (a, (3) (see Lemma\2^ below). 
According to J2E Theorem 2], this will imply thatj % ao w f3, i.e. {v\av2, V2f3v%, 
^17^3} is inconsistent. The consistency of the above set is, however, confirmed 
by our algorithm. Fig. [Tj] (left) gives a configuration of three regions {a, 6, c} 
which satisfies the constraints. 

We next show 7 is not decomposable w.r.t. (a, (3). First of all, it is easy to 
see that ct\ and ct2 in Eq. [39] are the only component single-tile relations of a 




in Eq. 38 
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ai= a 2 = 1 (39) 
\ 1 / \ / 

The next lemma characterizes when a basic relation is contained in the weak 
composition of a s and (3 for s = 1, 2. 
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Lemma 22. A basic relation tt is contained in a s o w f3 if and only if ir is a 
valid matrix whose component single-tile relations are also component single-tile 
relations ofn s , where tt s (s = 1,2) are matrices in Eq. 
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(40) 



Proof The proof is straightforward. We illustrate the weak composition in 
Fig. 15 (middle and right). Take ql\ in Eq. 39 as example. In the middle of 
Fig. 15 , there are six pixels marked "a." For each a, we have aa\b and bf3c. 
Moreover, dir(a, c) is a single-tile relation for each of the six a. These single-tile 
relations are components of 7Ti. The case of a 2 is similar and illustrated in the 
right of Fig. [l5j □ 
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Figure 15: Illustrations of evidences of weak compositions ao w /3 (left), a\ o w /3 
(middle), a 2 o w (3 (right) 

Then we arrive at our conclusion. 
Lemma 23. 7 is not decomposable with respect to (a, /3). 

Proof. Suppose 7 = 71 V 72, and j s C a s o w /3 (s = 1,2). This means, by 
Lemma [22} 7 S is a valid matrix whose component single-tile relations are all 7r s 's. 
For 1 < i, j < 3, we have = V j 2 [i, j], and j s [ij] < n s [ij} for s = 

1,2. Clearly, 71 [1, 1] = 71 [1,3] =0 and hence 72 [1, 1] = 72 [1, 3] = 1. Similarly, 
72 [1,2] = 72 [2, 2] = 72 [3, 2] = 0. Therefore, 72 is not 4-connected, hence not a 
direction relation matrix. This shows that 7 is not decomposable. □ 

The following example provides a set of three basic CDC constraints which 
is inconsistent but satisfies the condition of [32j Theorem 2]. 



Example 6. For basic relations a' \j3' ',7' as in Eq. 4J_ 

13' = 






(41) 
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we assert that j' is decomposable with respect to a' and f3' . Let a' s (s = 1, • • • , b) 
(the middle matrices of Eqs. "J^jty be the component single-tiled relations of 



It is easy to see j' s C a' s o w f3' (Eqs. ~J^^6). Evidences are given in Fig. 16 



where a configuration aa' s b, b(3'c, aj f s c is given for each s. Note that 7' is the join 
°f Is f or s = 1? " * " •> 5- According to Theorem 2 of \3$$ , 7' should be contained 
in a' o /3' ' , i.e. 

r' = {^1^2, v 2 f3'vz, ^17^3} 

should be consistent. Our algorithm, however, shows that V is not consistent. 
An argument is also given in Lemma\2J\ 
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(43) 
(44) 
(45) 
(46) 
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Figure 16: Illustrations of evidences of weak compositions a' s o w /3' (s = 1, • • • , 5) 



Lemma 24. T f = {via'v2, v 2 p f ^3, ^17^3} is inconsistent. 

Proof. Suppose a = {a, 6, c} is a solution to T f . We construct the frame S(a) 
and the cell set C(a) of a (see Dfn. [9|. Since (6, c) is an instance of /?', we 
know (M(b),M(c)) is an instance of l x (0') = (0 U fi) (0 U fi), i.e. 

(-M(fo), A4(c)) is an instance of one of the four basic RA relations 

0(g)0,fi(g)0,0(g)fi,fi(g)fi. 
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Take <g> as example. By a'[2, 2] = a'[l, 3] = a'[2, 3] = a'[3, 3] = 0, we know 
a° nbij = for (ij) = (22), (13), (23), (33). Similarly, by 7 '[2,1] = 7' [1,3] = 
7' [2, 3] = Y[3,3] = 0, we know a° nc fj = for (ij) = (21), (13), (23), (33). 
Therefore, the impossible regions separate the plane into two disconnected 
pieces, one is above the impossible regions, the other below (see Fig. 17). Since 
7^2,2] = 1, i.e. a° D A4(c) =^ 0, we know a has an interior point P which 
belongs to M(c), hence, the upper piece. By a' [3, 2] = 1, i.e. a° D 632 / 0, 
we know a has an interior point Q which belongs to the lower piece. This is 
impossible since we assume a is a connected region. □ 
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II III 



M(b) 



\ \ \\[ ! 
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Figure 17: Illustration of Lemma 24 



These two examples show that the weak composition algorithm proposed in 
[32] is not always correct. 



7 Consistency Checking for Two Variants of CDC 

The CDC algebra introduced in Dfn. [4] requires regions to be connected. This 
calculus has two variants in the literature. One, as introduced in [33], deals with 
cardinal direction relations between possibly disconnected regions; the other, as 
originally proposed in [11 , deals with simple regions, i.e. connected regions 
that are topologically equivalent to closed disks. In this section, we show our 
consistency checking algorithm designed for connected regions can be adapted 
to cope with these two variants. 

7.1 Cardinal Directions between Possibly Disconnected 
Regions 

For two (possibly disconnected) regions a, 6, similar to Dfn.|4j we write dir(a, b) 
— [dij]i j =1 for the direction relation matrix of a to 6, where dij is 1 if a°P\bij 7^ 0, 
and otherwise. A 3 x 3 Boolean matrix M is valid if there exist two region 
a, b such that M — dir(a, 6). It is easy to see that all but the zero 3x3 
Boolean matrices are valid. Each of these matrices represents a basic direction 
relation between possibly disconnected regions. We call the Boolean algebra 
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generated by these JEPD relations the Cardinal Direction Calculus for possibly 
disconnected regions, denoted as CDC^. 

Consistency checking in CDC^ is similar to that in CDC. Suppose M = 
{vi5ijVj}fj =1 is a network of basic CDC^ constraints. Similar definitions of 
regular solutions, meet-free solutions, and canonical solutions can be defined in 
CDCrf. Moreover, suppose a is a solution to J\f. We can transform a into a 
canonical solution a' = {a-}™ =1 of M. That is, a' is a regular, meet-free, and 
digital solution, and {Ix(a'i)}i=i and {I y (cLi)}2=i are canonical sets of intervals. 
Write rrii = I x (a>i) x Iy( a 'i)' We next show M has a maximal canonical solution. 
Actually, for each z, set Di to be the set of pixels in rrii that are disallowed 



for violating some constraint in M (see Eq. 24), and set bi to be the region 
obtained by deleting all disallowed pixels from rrii (see Eq. 25). We assert that 
b = {bi}^ =1 is the maximal canonical solution of N ' . 

Theorem 6. Suppose M = {viSijVj}fj =1 is a network of basic CDCd con- 
straints. If M is consistent, then b = {bi}^ =1 is the maximal canonical solution 
of J\f, where b{ is defined as in Eq. 



Proof. The proof is similar to that for Lemma 18 □ 



This shows, to construct the maximal canonical solution of a network of 
basic CDCd constraints, we need not to compute the connected components of 
bi. 

We next adapt our main algorithm to determine the consistency of a net- 
work J\f = {vi5ijVj}fj =1 of basic CDC^ constraints. As in the case of connected 
regions, Step 1 computes the projective IA networks J\f x and M y . If either is in- 
consistent, then J\f is inconsistent. Otherwise, Step 2 constructs their canonical 
interval solutions, and Step 3 computes bi for each i. In case M(bi) ^ rrii for 
some i then J\f is inconsistent. Otherwise, we go to the next step. The above 
procedures, as in the connected case, need at most cubic time. 

Since regions in a solution to CDC^ constraints are allowed to be discon- 
nected, we need not compute the connected components of each bi. Therefore, 
we go directly to Step 5, where we need check if dir(6^ , bj) = Sij holds for each 
pair of z, j. Suppose 5^ is represented as a 9-tuple (dfj)^ =1 . Similar to the 
connected case, we need only check whether b° D is nonempty for each <p 
with df- = 1. 

To check if b° D is nonempty, it is sufficient to show that there is a pixel 
p contained in bi D . Since bi may be disconnected, checking all boundary 
pixels is, however, insufficient. Note that checking all pixels in rrii H alone 
needs 0(n 2 ) time in the worst case for each pair z, j. This is undesirable since 
it will cost 0(n 4 ) time in total. 

We next show this could also be simplified. For a digital region b contained 



in T, write B(b) for the Boolean matrix that represents b (see Eq. 29), and 
define an n x x n y integer matrix M(b) as 

M(b) [k,l} = J2{B(b)[P,Q}-0<P<k,0<q<l} (47) 



36 



for each < k < n x and each < I < n y . It is easily to see that M(fo)[fc, I] is the 
number of pixels of b which are also contained in the rectangle [0, fc + 1] x [0, Z + l]. 

Given the Boolean matrix B(b), M(b) can be computed in 0{n 2 ) time by 
iteratively adding the fc-th column to the (fc + l)-th, and then iteratively adding 
the p-Xh row to the (p + l)-th. 

Lemma 25. Given a rectangle r = [x~ , x + ] x [y~ , y] and a digital region b, both 
contained in T = [0, n x ] x [0, n y ], b contains a pixel in r if and only if 
M(b)[x~ -l,y~ -1]+M(b)[x+ -l,y+ -1]> 

M(b)[x+ - 1,2T - 1]+M(b)[x- -l,y+-l], (48) 
where M(r)[-1,/] = M(r)[fc,-1] = 0. 

Proof Because M(b)[fc, I] denotes the number of pixels in b which are also con- 
tained in [0, fc + 1] x [0, 1 + 1], the number of pixels in b which are contained in 
the rectangle [x~ , x + ] x [y~ , y+] is M(b) [x + -l,y + -l]- (M(b) [x~-l,y + -i\ + 
M(b)[x + -l,y~- i\) + M(b)[x~ -l,y~-i\. The conclusion follows directly. □ 

Suppose T D = x [y~ ,y + ]. Since b{ C T, we have 

6* nmf = 6. n (rnmj) = b { n x [?r,2/ + ]. 
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By the above lemma, we know 6^ D mj contains a pixel if and only if Eq. 
holds, which can be checked in constant time. So, given M(^), the constraint 
Sij can be checked in constant time. Computing all M{bi) needs 0(n 3 ) time. 
Since there are n 2 constraints in total, Step 5, and thereby the who algorithm, 
can be finished in 0(n 3 ) time. 

Remark 2. For a consistent CDC^ network A/", the 0(n 5 ) algorithm proposed in 
[33] outputs a solution of M using possibly disconnected regions. It is difficult 
to extend this method to cope with cardinal directional constraints between 
connected regions. Actually, Navarrete et al. [23] tried to adapt this algorithm 
to cope with connected regions, and proposed an 0(n 4 ) algorithm for checking 
the consistency of a basic CDC network. The algorithm was based on Theorem 
1 in [23], which actually implied the following proposition. 

For a basic CDC network Af ', if all subnetworks of M 
involving four variables are consistent, then M is consistent. 

However, the counterexample constructed in Section 3.4 of this paper shows 
this is not true. 



7.2 Cardinal Directions between Simple Rgions 

In the definition of direction relation matrix (Dfn.|4|, we only assume connected 
regions. It is possible that these connected regions have holes. In the original 
work of Goyal and Egenhofer [H[|T2], objects are represented as simple regions, 
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i.e. regions that are tope-logically equivalent to closed disks. Interestingly, each 
direction relation matrix between connected regions can be realized by a pair of 
simple regions. Therefore, the set of cardinal direction relations between simple 
regions is the same as that between connected regions. We write CDC S for the 
qualitative calculus generated by these cardinal direction relations on the set of 
simple regions. 

In this subsection, we show this difference between simple regions and con- 
nected regions does not affect the consistency of cardinal direction constraints. 
In particular, we prove that a consistent network of CDC constraints always has 
a solution only using simple regions. The idea is to transform each connected 
region in the maximal canonical solution into a simple region without changing 
the relations between these regions. 

Suppose a is a digital region that is connected. Clearly, a has at most finite 
holes, where a hole of a is the closure of a bounded connected component of the 
exterior of a. It is easy to see that for a digital region all holes are digital and 
simple regions. 

Definition 17 (contact points). Let a be a digital region. A point P = (fc, I) 
is called a contact point of a if 

Pk,i Pk-1,1-1 C a <^ Pk-1,1 Pk,i-i 2 a - (49) 

Clearly, among the four pixels around a contact point, only two belong to a, 
and these two pixels are 8-neighbors, i.e. they are diagonally adjacent. 
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Figure 18: A digital region a with two contact points P = (5, 5) and Q — 
(5,6), where h\ and ft 2 are two holes of a, ftoo is the closure of the unbounded 
component of the exterior of a 



Suppose a = {a^}^ =1 is the maximal canonical solution of a basic CDC 
network. We next transform each ai into a simple region by deleting subpixels 



from a. Examples are show in Fig. 19 and Fig. 20 
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Figure 19: Transform a connected digital region without contact points into a 
simple region 




Figure 20: Transform a connected digital region with one contact points into a 
simple region 



Lemma 26. Suppose a = {a^}™ =1 is the maximal canonical solution of a basic 
CDC network. Let T be the frame of a, and let rrij be the mbr of aj for each j . 
If Pki is a pixel which is contained in a hole of ai, then there exists j ^ i such 
that pki C m j and a° D rrij = . 

Proof If p^ is contained in a hole if a^, then there exist pixels contained in ai 
that are, respectively, above, below, right of, and left oipki- That is, there exist 
ki < k < &2 and l\ < I < I2 such that Pk 1 hPk 2 i an d Pkh->Pki 2 are a ^ contained 



in ai. By the construction of a^, we know pki is a pixel in D{ (see Eq 24). This 



means, there exist j ^ i and (j) such that a° HraJ = and p^ C raj. If (j) ^ 5, it 
is easy to see ra^ also contains at least one of the four 'guarding' pixels Pk ± hPk 2 i 
and PkinPki 2 - This contradicts the fact that these pixels are contained in a^. 
As a consequence, we know = 5 and ra^ = rrr- = rrij. This means, there is 
j ^ i such that pki C rrij and a° D rrij = 0. □ 

Lemma 27. For two digital regions a, b, if a' C a and b f C b satisfy Eq. [5#] 
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and, respectively, Eq. 51, then dir(a,b) = dir(a',b'). 

M(a) = M(a f ) and for all pixels p, p C a implies (a') Dp ^ (50) 

M(b) = M.(b') and for all pixels p, p C b implies (b')° Dp ^ 0. (51) 

Proof This follows directly from the definition of direction relation matrix. 

Note that for any (j) = 1, • • • , 9 we have a D = if and only if a 7 D = 0, 

where = b^ = (b'Y . □ 



Recall a maximal canonical solution is a meet-free solution (cf. Dfn. 11). 
This implies that the mbrs of any two regions in a maximal canonical solution do 
not meet at a point. The following lemma is a consequence of this observation. 

Lemma 28. Suppose a = {ai}™ =1 is the maximal canonical solution of a basic 
CDC network. If P — (fc, V) is a contact point of a\, then among the four pixels 
around P, viz. Pk-i,i-i, Pkh Pk-i,i, Pk,i-i , one is contained in a hole of a±, two 
are contained in a\, the other is contained in neither a\ nor its holes. 

Proof. Without loss of generality, suppose pki and Pk-i,i-i are contained in 
a\. If Pk-i,i and Pk,i-i are contained in neither a\ nor its holes, then they 
are connected to a same pixel in the unbounded connected component of the 
exterior of a\. This means p^i is separated from Pk-i,i-i, which contradicts the 
assumption that a\ is a connected region. On the other hand, if each of Pk-i,i 



and Pk,i-i is contained in a hole of ai, then by Lemma 26 , Pk-i,i ^ TUj and 
Pk,l-i ^ mj' fo r some with a\ D rrij = and a\ D m f j = 0. Because p^ 
and Pk-i,i-i are not contained in rrij and m'^ we know rrij and m'j must meet 
at point P. This is impossible for maximal canonical solutions. Therefore, only 
one of Pk,i-i and Pk-i,i is contained in a hole of a±. □ 

By the above lemma, there are four different types of contact points. See 
Fig. [3j For convenience, we denote each type as a 4-tuple of symbols taken 
from {ft, a, For the four pixels, we start from the top left corner and go 
clockwise. These four types are written, respectively, as (ftaxa), (aftax), (xafta), 
and (axah). 



ft 


a 


a 


ft 


X 


a 


a 


X 


a 


X 


X 


a 


a 


ft 


ft 


a 



Table 3: Four types of contact points, where ft, a, x denote, respectively, a pixel 
that is contained in a hole of a, in a, and in neither. 



A contact point can be removed by deleting a sub-pixel from the a-pixel 
which follows the ft-pixel according to the sequence of the type of the contact 
point. Once all contact points are removed from each a^, the remaining holes 
are quite easy to cope with. 

Theorem 7. Let AT = {^iSijVj}^j =1 be a network of basic CDC constraints. If 
M is consistent, then it has a solution a = {a^}^ =1 such that each ai is a simple 
region. 
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Proof. Suppose a = {ai}^ =1 is the maximal canonical solution of AT. Then each 
cti is a connected digital region, which may have holes. We assert that, for each 
i, there exists a simple region a- such that a- C and .M(a-) = M(di), and 



dir(a^,aj) = dir(aj,aj) and dir(a 



dir(aj, aj) for any j ^ i. 



To prove this statement, we first subdivide each pixel into 25 equal sub-pixels 



(see Fig. 21). 
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Figure 21: Subdivision of a pixel into 25 sub-pixels 



For a contact point P = /) of a^, without loss of generality, assume P has 
type (haxa). Removing the 1/25 sub-pixel which contains P from a^, we obtain 
(after necessary regularization^]) a new region (see Fig 22). This procedure can 
be applied to all contact points of at the same time. Write a* for the resulted 
region. Since each pixel p of has at most four contact points, the revised 
region a* contains at least 21 of the 25 sub-pixels of p. This implies that 
a* C en and M(a*) = M(di). It is routine to check that a* is still connected 



but has fewer contact points as well as fewer holes. By Lemma 27 we know 
dir(a*,<2j) = dir(a^,aj) and dir(<2j,a*) = dir(aj,a^) for any j ^ i. 
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Figure 22: Remove a type (/iaxa)-contact point 



Now, we have removed all contact points of a^. This means, any two holes of 
a* are disjoint, and each hole of a* is disjoint from the closure of the unbounded 
component of the exterior of a*. For each hole h of a*, select a pixel p{h) in 
h which has the highest y-mdex but whose left 4-neighbor is out of h. We cut 



1 Note here by regularization we mean the topological regularization of a set a, which is 
defined as a° . 
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a slot from a* to connect the hole h with the exterior of a*. For a pixel p 
above p(h), if p and all pixels between p and p(h) are contained in a*, then 



delete the middle column of sub-pixels from p (see Fig. 19). After necessary 
regularization, we obtain another connected region which has fewer holes than 
a*. Applying this operation to holes of a* one by one, we obtain a connected 
region aj, which has no holes. That is, a[ is a simple region. By Lemma 27, we 
know dir(a-,<2j) = dir(a^,aj) and dir(<2j,a-) = dir(aj,a^) for any j ^ i. 

We claim {a-}™ =1 is a solution of M which consists of simple regions. For 
any i ^ j, because M(a'j) = M(aj) and dir(a-,<2j) = dir(a^,a J ), we know 
dir(aj,a^) = dir(a^aj). The conclusion holds already. □ 

As a conclusion, we know the satisfaction problem over CDC S can be decided 
in the same way as that over CDC. 

Theorem 8. The consistency of a network of basic CDC S constraints can be 
determined in cubic time. 



8 Conclusion 

This paper provided a cubic algorithm for checking consistency of basic CDC 
networks, which was earlier observed as impossible [23 for connected regions. 
If a basic CDC network is consistent, our algorithm also generates the maxi- 
mal canonical solution. This general algorithm was then applied to solve two 
subproblems: the pairwise consistency problem and the weak composition prob- 
lem. Through careful examination, we showed that the cardinal direction rela- 
tion model used in [3] is not coherent, and showed by examples that the weak 
composition algorithm obtained in [32 is not always correct. 

Although devised to solve cardinal directional constraints between connected 
regions, our main algorithm can also be adapted to cope with cardinal directional 
constraints between possibly disconnected regions as well as those between sim- 
ple regions. For a basic network of CDC^ constraints, our algorithm determines 
in cubic time if it is consistent. Compared with the 0(n 5 ) algorithm reported in 
[33] , our algorithm is more efficient. As for cardinal direction constraints over 
simple regions, we proved that a basic CDC network has a solution using only 
simple regions if it is consistent. This suggests that CDC does not distinguish 
between simply connectedness and connectedness. 

Most potential applications of qualitative spatial reasoning require multiple 
aspects of space. Combining spatial constraints of different calculi is a very im- 
portant problem in the research of qualitative spatial reasoning. Some work has 
been done in this direction (see e.g. [9j [16j [22] ) . In particular, [22 points out 
that reasoning with basic RCC8 and basic RA constraints is in P, but reasoning 
with basic RCC8 and basic CDC^ constraints is NP-Complete. Note that in 
RCC8, RA and CDC^, spatial variables are interpreted over possibly discon- 
nected regions. It is still open if the above results for possibly disconnected 
regions still hold for connected regions. In particular, we do not know if rea- 
soning with basic RCC8 and basic CDC constraints is still decidable if spatial 
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variables are interpreted over connected regions. It is worth noting that con- 
nectedness is a tough requirement as far as topological relations are concerned. 
Although path-consistency suffices to determine if a basic RCC8 network is 
consistent, it was recently proved that deciding if a basic RCC8 network has a 
solution using connected regions is NP- Complete [30] . 
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Appendix A: Proof of Lemma [6] 

Lemmapl Suppose S = [dij]^j =1 is a basic CDC relation. Then the x-projective 
interval relation l x (S) is the IA relation associated to the vector (^1,^2,^3); i.e. 
(7, J) G t x (S) iff dir(I,J) = (di,d 2 ,d 3 ) ; where dj is if ' Ef =1 dij = and 1 
otherwise. 

Proof. We need to prove two things. First, suppose (a, b) G S. We need show 
v(/ IC (a), I x ip))i the direction relation vector of (I x (a), I x (b)) is (^1,^2,^3). Sec- 
ond, suppose v(/i,/2) = (^1,^2^3)- We need show there exist two connected 
regions a, b such that (a, b) G 5 and I\ = I x {a), I2 = I x (b). 

The first part can be proved as follows. Set I x (a) = [x~(a), x + (a)], I x (b) = 
[x~(b), By the definition of cardinal direction matrix, x~(a) > x~(b) 

if and only if a° D bn = 0, i.e. Sn = 0, for i = 1,2,3; (x~(a), x + (a)) D 
(x~(b), x+(b)) = if and only if a°n^2 = 0? ^2 = 0, for z = 1, 2, 3; x + (a) < 
if and only if a° D ^3 = 0, i.e. ^3 = 0, for i = 1,2,3. Therefore, the 
vector (di, ^2,^3) as defined here is the direction relation vector of (I x (a), I x (b)). 
In other words, (I x (a),I x (b)) is an instance of t x (S). 

To prove the second part, we note that 5 as a direction relation matrix is 
4-connected (see Prop.[2|. The proof is by construction. Take 




as example. In this case (^1,^2,^3) = (1, 1, 0). Suppose dir(/i, I 2 ) — (di, g?2? ^3)- 
Set I\ = [u, v], I2 = [s, t]. Then we have u < s < v < t. There are two subcases. 
If v < t, put ii and /2 on the x-axis of the orthogonal basis of the plane. Set 



a, b to be the polygons in Fig. 23 (left). It is clear (a, b) G S hold. The case of 
v = t is similar and illustrated in Fig. [23] (right). The same method applies to 
all the other basic CDC relations. □ 
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Appendix B: Proof of Lemma [9] 
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Figure 23: Construction of Cardinal Direction Relations from x-projective rela- 
tions 



Lemma pi An IA network M = {viHjVj}fj = i of constraints over 6* 



int 



consistent if and only if M = {viiijVj}™j = i is satisfiable. 



Proof We first observe that AT is a refinement of J\f. Therefore, if Af is consis- 

~]}it=i is a solution to 



\u 4 ,11, 



tent, so is AT. On the other hand, suppose {Ii 
N ' . We construct a solution to J\f . 

We call a point u G M = {u~ ,uf}f =1 a meet point if there exist i ^ j 
such that u~l = uj , i.e. Ii meets Ij. We use induction on the number K of 
meet points. If there is no meet point, {Ii = [u^ , uf]}f =1 is also a solution 
to M. Suppose the result holds for K = m. We show this also holds for 
K = m + 1. Suppose w = u~ Q is the largest meet point. We define a new 
solution {Ji = [Vi~ , v^~)}2 =1 to M which has fewer meet points. Set 



and 



w - 



■e/4, 



if u- = w, 
otherwise. 



where e is the smallest distance between different points in M. We next show 
each pair (Ji,Jk) is also an instance of i^. To this end, we need only consider 
the case when either Ji ^ Ii or J^ ^ Ik- Note that Ji ^ Ii if and only if u~[ = w. 
This means we need only consider the cases when u\ 



u k , and v i 



Case 1. If u i = w and u k = w, then v. 



e ,w + 6;/4}, v k e {u k , 
Vu = w 



= w or u k = w. Recall 
w + e/4}. 

+ e/4. The ordering of 



is the same as that of u- , u] 



Case 2. If u i = w and u k = 



,u k: u k . 



w, then Uu < u 



+ _ 



u- < vj and v, < vt < 



Case 3. 

JiPJk- 

Case 4. 
as that of 



~. This means Ik^h and JkPJi- 



If 



If 



w and uj~ = w, then, similar to Case 2, we have iiini/c and 
= w and u k ^ w, the ordering of v~,v^,v k ,v k is the same 
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Case 5. If = w and uf ^ w, the ordering of v~,v^,v^,v^ is the same 
as that of 

So {Ji}f =1 is also a solution to A/". Note that either Case 2 or Case 3 is 
true, because w is a meet point. Moreover, w is not a meet point in the new 
solution anymore since all 'meets' instances at this point has been changed to 
'before' instance. Therefore, {Ji}f =1 has fewer meet points than {ii}f =1 . By 
the induction hypothesis, we know M has a solution. This ends the proof. □ 
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